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SUMMARY 


A design criterion for input functions in laboratory tracking 
tasks resulting in efficient parameter estimation is formulated. 
The criterion is that the statistical correlations between pairs 
of parameters be reduced in order to minimize the problem of non- 
uniqueness in the extraction process. The effectiveness of the 
method. is demonstrated for a lower order dynamic system. 


INTRODUCTION 


The art of human operator modeling has progressed consider- 
ably since the pioneering work of Sheridan (ref. l), Elkind et . 
al . (ref. 2), and McRuer et . al . (ref. 3). Many of the accepted 
transfer functions for the human operator in compensatory track- 
ing tasks are given in reference U. The model structure for the 
human operator is not known as precisely for such systems as 
the aircraft, for example. There are large variations from one 
subject to another and from run to run. In most cases, the pilot 
model includes such physical param’eters as 'a static gain, an 
effective time delay made up of transport delays, and high frequency 
neuromuscular lags or leads and turns which represent the low fre- 
quency characteristics of the neuromuscular system dynamics. In 
addition, a remnant function is generally included in the model 
to account for pilot anomalies and unsteady behavior. It has been 
pointed out by various authors (refs. 5-6) that large remnants and 
the quasi-predict able nature of the inputs cause difficulties in 
the extraction of the physical parameters in the pilot model. 

Early attempts to determine the best set of parameters, in particu- 
lar pilot models utilizing pilot response data, relied upon analo’g 
matching techniques (ref. 7). More recent techniques have applied 
Kalman filtering methods (ref. 8), maximum likelihood methods 
(ref. 9)j or Newt on-Raphson methods (ref. 9) to the system identifi- 
cation p r Ob 1 em . 
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One of the practical difficulties in the extraction of system 
parameters from experimental data is the non-uniqueness in the 
solution for the parameter values (ref. 10 ) . In the case of 
human operator dynamics, the non-uniqueness problem is amplified 
by the high correlations which can exist between the pilot's effec' 
tive time delay and the time lags associated with the equalization 
characteristics' discussed by McRuer et. al . (ref. 3). Rormall^, 
the time delay is either assumed known or represented by a Fade 
form which places it in the role of either lead or Tag constants 
( ref . 9 ) . 

This paper addresses the task of reducing the non-uniqueness 
and possibly the effect of the remnant by proper design of the 
forcing function or disturbance used in laboratory designed track- 
ing t asks . 


MODELS 

The performance of a human operator in many tracking tasks 
can be modeled adequately by a quasi-linear describing function, 
The describing function model consists of a transfer function 
Y (s) and a remnant n(t) as shown in figure 1 by the block di a- 

gram of a typical compensatory tracking task. 

One of the more generalized transfer functions for the com- 
pensatory control tasks is discussed in reference U and is writ- 
t en as 


Y (s) = K e 
P P 


or (T s+1) 
-ST L 


(TjS+1) 


(Tj^s + 1) 


{T^s + 1)(T„^.-M)[{^)% ^ s.ll 


( 1 ) 


In such a model, the parameters K , x, 

and are generally poorly determined 

using either analog matching techniques 
methods in conjunction with pilot respo 

Depending upon the plant to be con 
meters in equation (l) can be eliminate 
For purposes of illustration of how the 


T. 


^k' 


•Nl' 




L' "I' "k' 
and are improved upon by 
or parameter estimation 
nse data. 


trolled, some of the para- 
d from consideration, 
problem can be put into 
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a state space formulation for wliich the methods of modern estima- 
tion algorithms can he utilized, the following special case will 
he considered: 


i 

li 


Y 



Y 

P 


(s) 


K 

2 

s 


K^(T^s + l) e"®"^ 
(T^s+l) (T^s+l) 


( 2 ) 

(3) 


In addition, it is assumed that the remnant function is the re- 
sult of white noise w (t) filtered through a second-order linear 
filter according to ^ 


n(t) + n(t) + n(t) = w^(t) 
From the block diagram 


(U) 


u(t) + (T^+T 2 )u(t) + u{t) = K^[T^ e(t-t) + e(t-x)] (5) 

c (t ) = K6 (t ) (6 ) 

Equations (h)-(^) can he written in terms of the state variables 

T 

y(t) = (y^,y^ 


= [ c( t ) , c (t ) ,u(t ) , r(t) + u(t + x) + T u(t + x), 

1 

n (t ) ,n (t ) 

as 

y(t) = y(t) +F^ y(t-x) + h(t) + d(t) 


( 7 ) 

( 8 ) 
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where 
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The measurements used in the estimation procedure are gener- 
ally such quantities as the plant output c(t), the stick output 
6(t), etc. These quantities are written as a vector 

z(t) = h(y) + v(t) ( 13 ) 

where v(t) is assumed to he a white noise process with known 
statistics . 

DESIGN METHOD 

The design of optimal inputs for non-time delay differential 
systems has been investigated by numerous authors. Goodwin 
(ref. 11 ) designed an input which minimized the covariance of 
the error estimate. Mehra (ref. 12) designed an input to maxi- 
mize the sensitivity of the system output to the system para- 
meters. The philosophy adopted in the present work is to design 
an input to minimize directly the correlation between certain 
pairs of parameters. 

The problem is stated in terms of a di f f erent i al- di f f erenc e 
equation as given in equation (8) restated as 

y(t) = f(t, y(t), y(t-x), p) + u(t) + w(t) ( l4 ) 

with initial values 

y(t) = 0,t -T<t<t 
o o 

X is the constant time delay and p is the unknown parameter 
vector. The vector u(t) is to be chosen optimally based on 
a design criterion and w(t) is the white noise process related 
to the remnant . 

The system measurement from which p is estimated is taken as 
the system state vector y(t). The noise in the measurement is 
assumed a Gaussian white noise process with zero mean and co- 
variance R. 

In addition to the state equation expressed by equation (lU) 
it will prove useful to introduce additional state variables 
defined by the elements of the sensitivity matrix A(t) and the 
time dependent covariance matrix of the error in the estimate 
C(t) defined respectively according to 

A(t) = 8 f A(t) + 9f A(t-x) + ^ ( 15 ) 

3x ( t ) 9x ( t -X ) 9p 
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A(t) = Of t -T<t<t 

O _ _ O 


and 


C(t) = -C(t) A^(t) R"^ A(t) C(t) 


(16) 


C(t^) = A^(t) R“^ A(t)dt] 


-1 


(IT) 


With these additional ec[uations , it is possible to incorporate 
equations (l4)-(l6) into one differential equation of the form 


F(t, X(t),X(t-T), X(t ),x(t-T),u(t) , w(t), p) = 0 
where the augmented state vector X is defined as 

X(t) = [x^(t) 
and 
Ai) 


... 




(t) — [C(t)^ C(t)/.... C (b)] 
-LX d.d. mm 


T 


c<^>(t) = [c^3(t), 

« • • • 

C^”^t) = c^^(t) 


(18) 


.A^^^^(t)]'^ 


( I 1 3x (t ) , 9x (t ) , . . 

A^^At) = [ i ■— i- 

3Pl 3?2 


3x (t ) „ 

3p ^ 
m 


/ ) 3x (t) 3x (t) .... 3x (t) 

^^1 ^^2 3Pm 

The optimal input is designed to minimize the average value 

of a performance index representing a weighted sum of the squares 

of the correlation coefficient plus an arbitrary function in the 

control variable* The correlation coefficient between parameters 

p. and is defined as 

1 J 


The performance index can then be written in conventional 
form as 
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t - 

J(u) = E {Mx(t^)) + L(x(t), u(t),t)dt} (19) 

where the first term represents the weighted sum of the squares 
of the correlation coefficients . The second term can he chosen 
according to whatever physical requirement is desired between the 
state and control vectors. 


EXAMPLE WITH LOW ORDER DYNAMICS 


As an example of the design concept, consider the scalar 
system 

x(t)+^x(t-T) = u(t) + w(t)/0<t<^+ (20) 

x(t) = O, -T<t <0 

where 5 and x are the system parameters to be estimated and u and 
w are the input and white noise funtions respectively. The in- 
put or control is constrained according to 

|u| < 2 

and is assumed zero prior to time zero. 


The performance index, which is the square of the correlation 

T is 


coefficient between parameters 5 and 

ai (t ) a2 (t )dt ) 


J ( u) = 


(/o 


(t)dt) (/^ a^‘=^(t)dt) 


( 21 ) 


The form of the optimal control is a bang-bang control expressed 
as 

u(t) = 2sgn g(x(t,x),t) ( 22 ) 

where g is the switching function which can be computed from the 
maximum principle. The optimal estimate of x(t-x) denoted 

x(t,x) has been shown by Kwakernaak (ref. 13) to be of the form 

x(t,x)=/^K°(t,x,a)z(a)da (23) 

where the Kernel function K° is a function of the state equation 
and noise source . 

To continue, a further simplification to the deterministic 
ca?ie will be made. this case, the sensitivity matrix is the 

two component row vector 
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( 24 ) 


(25) 


By use of Laplace transform techniques, the solutions to 


equations (20) and ( 25 ) are written as 

x(t) = h(t-cr) u(a)da ( 26 ) 

ai(t) = -x(t) + 1 /q h(t-a)x(a)da ( 27 ) 

a 2 ("t) = h(o)x(t-T) + h ( t-ox ( a-T ) do (28) 

where h(t) is the unit impulse response 

h(t> = I e ^ [1-T/^e (29) 

and X is the kth eight value of the ch ar ac t eri st i c ’ e quati on 
^ _ Xt 

+ e = o ( 30 ) 


The effectiveness of the design technique for ^ = t = ll2 
is shown in figure 2. The correlation is reduced from a value 

of 0.60 to 0.09 Ly means of the design technique. 

CONCLUDING REMARKS 

An optimal design technique for disturbance functions used in 
laboratory tracking tasks has been developed. The objective is 
to reduce statistical correlations between various parameter 
pairs to effect a more efficient parameter extraction from noisy 
experimental data. Preliminary calculations indicate the design 
methods to be effective. 


A(t) = r 3 x( t ) # 9 x( t ) 1 

■•95 9a ■' 

which satisfy the equations 

ai(t) + lai(t-T) = -^x(t) 

K 5 

a . 2 ( t ) = l^x ( t - T ) 

5 
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